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Nomenclature 


A area, 

abs absorbed 

C specific heat at constant pressure, J kg~^ 

Dj inner diameter of pipe, m 

Do outer diameter of pipe, m 

h heat transfer coefficient, W m~^ 

ij,k indicator of i, j and k^^spatial position in x,y and z directions respectively 

k thermal conductivity, W m”^ 

It solar radiation flux on the tilted collector per unit area, W m“^ 

m total fluid mass flow rate, kg s“^ 

n number of nodal points 

q heat flux, W m“^ 

Q rate of heat transfer, W 

P perimeter, m 
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S amount of heat flux absorbed by the absorber plate, W 
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V volume, m^ 

F’ collector efficiency factor 
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UTfi,UTja uncertainty in measurement of inlet and outlet temperature 

UvjUt uncertainty in measurement of volume and time 

Vm uncertainty in mass flow rate 

U,, uncertainty in efficiency 

Bj biot number 

x.y,z spatial coordinates 


Greek Symbols 

a absorptivity of absorber plate 

0 collector tilt, radian 

7 surface azimuthal angle, radian 


5 


thickness, m 



NOMENCLATURE 


vii 

emissivity 
latitude, radian 
reflectivity; density, kg m“^ 

angle of incidence of radiation on a tilted surface, radian 
refraction angle for cover2, radian 
refraction angle for cover 1, radian 

angle of incidence of radiation on a horizontal surface, radian 
transmissivity of glass cover 
Stefan-Boltzmann constant, 5. 67*10“® W m“^ K“'* 
efficiency 

1 * 10-5 

Subscripts 

ambient 

beam 

cover2 

coverl 

diffuse 


fluid 



menclature 


viii 


insulation, inlet conditions 
mean 

output conditions 
absorber plate 


wind 



Abstract 


The objectives of present work were to develop a multidimensional transient model 
of a flat plate solar collector to predict its performance and conduct experiments to 
validate the model. 

The performance of flat plate solar collector depends upon solax insolation, latitude, 
time of day, angle of tilt, mass flow rate and inlet fluid temperature. The model devel- 
oped here can incorporate the effect of all the above factors to predict the performance 
of a flat plate solar collector. The flat plate solar collector was designed, fabricated and 
tested on different days with different mass flow rates, angles of tilt and with varied 
distance between the absorber plate and the cover glass plate and also that between 
the two covers in a double cover design. The model was validated with the experiments 
conducted on the collector. 

The model equations were developed considering different elements lihe covers, ab- 
sorber plate, insulation, fluid etc. which were coupled in a global algorithm. These 
equations were solved using finite difference technique and computed results showed 
that (a) Rise in fluid temperature was maximum between 12 noon and 1P.M. (b) PMte 
mean temperature was maximum between 1 P.M. to 2 P.M. (c) The mean efficiency 
of collector decreased as the inlet fluid temperature was increased, (d) The mean effi- 
ciency was constant for flow rate greater than 40 kg/h which indicated the minimum 
flow rate for the collector, (e) The theoretical results were in good agreement with the 
experimental data from 9 A.M. to 2 P.M. which implied that transient model developed 
here was validated by the experimental results. This model therefore could be used to 
characterize a flat plate solar collector as the function of all the relevant parameters. 



Chapter 1 


INTRODUCTION 


1.1 Need for alternative source of energy 

The world population is growing fast. People around the world are demanding more and 
more energy as they aspire to raise their standard of living which is directly associated 
with the energy consumed. Nations consuming more energy on a per capita basis 
have better living stmidard than the others. Because of this natural desire to improve 
living standards, world energy consumption is increasing at a faster rate than world 
population [ 35]. 

Energy is the major infrastructure requirement for any developing or developed 
country. Today most of the world’s energy demand is met by fossil fuels, mainly 
petroleum and natural gas. Enormously growing demand for energy and increasing 
exploitation of the available conventional energy resources is causing a rapid depletion 
in their reserves, which may thus not last for a long time. Conventional energy sources 
are polluting and leading in global pollution and global warming. There has been 
growing interest in the development of alternative energy sources because of global 
pollution, rising prices and limited reserves of fossil fuels. Solar energy, wind energy, 
biomass energy and tidal energy etc. are few of the available alternative energy sources. 
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While the potential for energy generation from these sources is immense, one of the 
promising option is to make more extensive use of energy derived from the sun. 


1.2 Solar energy 

Taking a global energy perspective, it is felt that solar energy is one of the most 
promising alternate sources because of its world wide availability particularly in coun- 
tries lying within 35° N and 35° S latitudes. Solar energy has a great future in the 
tropical country like India, where many places receive a peak solar flux of the order of 
1 kWfm^ [19). In principle, solar energy could supply all the present and future energy 
needs of the world on a continuing basis. Solar energy, unlike fossil fuels and nuclear 
power, is an environmentally clean source of energy. Another advantage is that it is 
free of cost and available in abundance in a large number of countries of the world. 

There are also problems associated with proper utilization of solar energy. It is 
intermittent, low density and conversion to the electrical power is not very efficient. 
Solar energy can be used both directly and indirectly. It can be used directly in a 
variety of thermal applications like heating water or air, drying, solar distillation and 
cooking. The heated fluids can in turn be used for applications like power generation 
or refrigeration. A second way in which solar energy can be used directly is through 
the photovoltaic effect in which it is converted to electrical energy. Indirectly, the 
sun causes winds to blow, plants to grow, rain to fall and temperature differences to 
occur from the surface to the bottom of the oceans. Useful energy can be obtained for 
commercial and noncommercial purposes through all these renewable sources. 
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1.3 Solar collectors and their main components 

In the last few decades, new solar collectors of different shapes have come up. The 
three most prevalent collector types are fiat plate, evacuated tube, and concentrating 
type. The flat plate collectors have evoked considerable interest because they are the 
simplest among the solar energy collecting devices, least expensive and have a wide 
range of potential for heating and cooling applications. A schematic diagram of a 
fiat plate solar collector is shown in Fig. 1.1. They usually consists of the following 
components : 

1. Absorber plate 

An absorber plate is generally a metal sheet of high thermal conductivity. It is coated 
with black paint or given a special coating so that it absorbs the incident solar radiation 
efficiently and loss of heat by radiation from the absorber plate is minimized. An 
absorber plate absorbs the incident solar radiation, converts it into heat and conducts 
the heat to the fluid passages. 

2. Cover plate 

A gla.ss sheet of good quality which is transparent to the incoming solar radiation is 
fixed above the absorber plate. This prevents convective heat loss from the absorber 
plate and prevents infra-red radiation from the plate to escape to the atmosphere. The 
cover plate may be one or more in number. 

3. Fluid 

1 he heat transfer fluid is required to remove the heat from the point of collection and 
transferring it to the point of application. It is a means to extract the heat absorbed 
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by the absorber plate. This hot fluid may directly be used for heating purposes or may 
transfer its heat to another secondary fluid through a heat exchanger. 

4. Insulation 

The bottom and side part of the absorber plate is separated from the ambient conditions 
by means of a layer of insulation materiaJ. This insulation prevents the loss of heat 
by the plate to the surroundings. Usually the insulation used is fiberglass, glass wool, 
wood shavings, Styrofoam, saw dust, urethane, therraocole etc. 

5. Collector box 

An enclosing box whose principal functions is to hold the other components of the 
collector in their relative positions and to protect the collector plate and insulation 
material from the weather. 
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1.4 Scope of present work 

The objectives of the present work are ; 

• 'ib (lovelop a imilticlimenHlonal transient model of a Hat plate solar collector to 
predict its performance. 

• To conduct experiments to validate the model. 

The performance of a flat plate solar collector d^ends upon the solar insolation, lati- 
tude, time of day, angle of tilt, mass flow rate and inlet fluid temperature. The model 
developed here attempts to incorporate the effect of all these factors and is validated 
by experimental data. 
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LITERATURE SURVEY 

2.1 General remarks 

The solar radiation reaching the earth’s surface in any location varies considerably 
with atmospheric conditions, hour of day and season of year, with the result that solar 
heat collectors always operate in an unsteady mode. The absorber plate acts as a fin 
to conduct heat to working fluid in collector tubes, implying the existence of thermal 
gradients in the plate near the tube, in direction normal to the fluid flow. Hence, at 
any time, the absorber plate also operates in a nonisothermal mode. 

The transient modeling of a collector in a solar heating system is important since 
an accurate calculation of the collector’s performance provides information for the 
assessment of the collector system for the systematic investigation of collector d^ign 
and for understanding of solar energy contribution in the passive heating of buildings 
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2.2 Review of published work 

The first noticeable work available in literature in connection with flat-plate collector 
was by Hottel and Woertz [14]. They developed one dimensional steady state model 
for thermal behavior of a solar collector. The interaction between the incoming solar 
radiation absorbed by the covers and outgoing loss was neglected. Gutierrez et al. [11] 
used hybrid computer to simulate solar water heating systems using specified storage. 
They examined the effect of three types of auxiliary heating systems with three tanks 
of different size at various time distribution and loads. Klein et al. [17] investigated 
the effect of thermal capacitance in the modeling of the performance of a flat-plate col- 
lector. They investigated the three different models, a quasi-steady-state model with 
zero capacitance, second one with a single value of thermal capacitance and a third 
model dividing collector into many isothermal segments, or nodes. They showed that 
for hourly meteorological data, the zero capacitance model was adequate. Klein et al. 
[18] described the simulation of dynamic behavior of transient solar energy system by 
solving sets of differential and algebraic equations. De Ron [8] developed the model for 
a flat-plate solar collector considering dynamic effect. His verification and estimation 
procedure for the model proved to be a useful tool for comparing various collectors. 
Morrison and Ranatunga [21] investigated the response of thermosyphon solar water 
heaters to step changes of insolation. They developed a mathematical model to simu- 
late the transient performance and measurement of the transient flow development in a 
thermosyphon circuit. Hobson and Norton [13] conducted detailed analysis for the heat 
transfer and fluid flow within direct thermosyphonic solar energy water heaters. The 
collector energy equations were cast in a two-dimensional form in order that heat trans- 
fer and thermal capacitance effects were simulated accurately at the small flow rates 
encountered commonly in these buoyancy driven systems. Saito et al. [28] presented 
numerical solutions for a transient one dimensional model pf a solar water heater. Pier- 
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son and Padet [26] described a method to calculate the time constant of solax collector. 
They proposed an analytical expression based on experimental observations and on a 
thermal analysis in unsteady state. They also studied the influence of collector param- 
eters on time constant. This helps in comparing the transient performances of different 
solar collectors. Chiou [5] presented a two-dimensional steady state model to study 
the effect of nonuniformity in the fluid flow distribution on the performance of a flat 
plate solar collector. Oliva and Costa [24] presented a numerical simulation of solar 
collectors in which they considered the effect of nonuniform and unsteady boundary 
conditions. They also studied the influence of shadows due to external obstacles and 
cover supports to the collector. Zhao et al. [36] presented a computer model derived 
from the basic laws governing thermal energy exchanges between surfaces for transient 
simulation of flat plate collectors. The approach of subdivision along the collector was 
used in the development of energy equations on each of the surfaces in the collector. 
The model was validated with experimental data. Onyegegbu and Morhenne [25] pre- 
sented an unsteady two-dimensional analysis of a flat-plate solar collector subjected 
to time varying insolation with considerable diffuse components as a better character- 
ization of practical solar collectors. Their ’analysis considered the thermal masses of 
absorber, tube, glazing , and the working fluid in the system and also the different 
optical and thermodynamic properties of beam and diffuse radiation. Srivastava and 
Kant [30] studied multidimensional transient analysis for single glazing . 

Iqbal [15] studied the effects of free convection superimposed on forced flow in 
inclined tubes of a solar collector. The theoretical analysis considered fully devel- 
oped Laminar flow under uniform heat flux condition and constant pressure gradient. 
Choudhury [6] discussed simple and efficient parameterization of insolation under par- 
tially cloudy sky and compared it with a set of exact radiative transfer results. The 
albedo of the earth-atmosphere system was also formulated and compared with a set 
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of exact radiative transfer results. Gaiii and Symons [10] studied the simulation of the 
influence of cover design on the thermal performance of flat-plate solar collectors for 
use at temperature of 150°C. They presented detailed results of the effect of changes 
in cover materials, cover surface treatments, cover system configuration and absorber 
plate surface treatment on both the instantaneous efficiency and the long term solar 
contribution of flat-plate collectors. They showed that flat plate collectors of single 
high transmittance cover, a convection suppressing device such as a honeycomb and 
selective absorber surface yielded long term solar contributions comparable to those of 
evacuated tubular collectors where as other simulated configurations (single or multiple 
cover systems) were significantly inferior. Bhide et al. [3] presented a simplified method 
for the comparison of thermal performance of two similarly fabricated flat plate col- 
lectors with absorber coated with coatings of known solar absorptivity and emissivity. 
They essentially determined the insolation values at which the two collectors gave the 
same thermal performance. Curves were presented which could be directly utilized by 
the designer for selecting the coating having better thermal performance. Thomas et al. 
[33] analyzed angular response measurement for four different types of collector, each 
type tested in three different laboratories. An analysis of the measurement procedure 
showed that experimentally determined angular response parameters were subjected 
to relatively large uncertainties. Morrison et al. [20] investigated the effect of load 
draw-off pattern and auxiliary boosting on the performance of vertical and horizontal 
tank solar water heaters. The primary factor governing the dependence of performance 
on load pattern in single tank systems was conduction between the auxiliary boosted 
zone and the solar preheat zone. A method for quantifying the effect of conduction in 
solar tanks was presented. Behnia and Morrison [2] investigated free convective flow 
and thermal structures in an open ended inclined cylindrical thermosyphon with the 
application to an evacuated tubular solar collector. The flow structure was visualized 
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using dye and rhescopic tracer particles and temperature measurements were made by 
a traversing thermocouple rake. Garg [12] studied the effect of dirt on the transmit- 
tance of solar radiation through various glass plates and plastic films, which are used 
as a transparent cover for flat-plate collectors. The dirt correction factor for glass plate 
inclined at an angle of 45 deg from horizontal was 0.92. The correction factor was 
greater for the plastic film as compared to that for the glass plate for any inclination. 

Russell and Guven [27] studied all-plastic liquid flat-plate collector. Such a collector 
eliminates the need for coatings, metals, and all other materials except plastics. The 
concept appeajs to offer possibilities for the mass production of a low-cost, lightweight 
solar collectors which will be reasonably efficient. Mathematical models were developed 
for study of this concept. A model panel was fabricated, and computer r^ults were 
compared with the experimental data for the efficiency of the panel. Nayak et al. [22] 
studied on solar concrete collectors used for providing domestic hot water. Although 
they are less efficient, yet such collectors integrated with building structures are likely to 
be more cost effective than conventional flat-plate collectors which use metallic absorber 
plates. The collectors axe made up of thin concrete slabs with a network of PVC tubes 
embedded inside concrete. A layer of galvanized iron wiremesh on either side of the 
PVC tubes provides the necessary reinforcement to the concrete. The top of the slab 
was painted black and glazed, while the back insulation was made of autoclave cellure 
concrete. Nordgaard and Beckman [23] presented a new method for predicting the 
radiative transfer in an absorbing and isotropically scattering, non gray plane-parallel 
atmosphere. The technique was based upon the F (“F-hat”) concept and developed 
to determine the radiative transfer at short wavelengths for the transparent insulation 
material monolithic silica aerogel. The performance of a system using aerogel collectors 
was compared with system using other types of high performing collectors. 

Jenkins and Hill [16] compared flat plate water-heating solar collectors using the 
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BSE and ASHARE test procedures. Five solar collectors were tested according these 
procedures and the results were compared. These collectors were modular, flat-plate, 
water-heating, and included single and double-glazed designs with or without selec- 
tively coated absorber plates. In both procedures, collector efficiency curves were 
plotted. The ASHARE procedures consisted exclusively of outdoor testing whereas 
the BSE procedures required a combinations of outdoor and indoor testing (no irradia- 
tion) to determine the collector’s optical and thermal loss characteristics, respectively. 
During the indoor testing in this study, the environmental test conditions were con- 
trolled and regulated by use of specially built environmental simulators and thus the 
effect of wind speed and “sky” temperature on the thermal loss characteristics of the 
collectors were investigated. 

Amar et al. [1] presented the review of the literature on the test methods of flat-plate 
solar collectors under varying environmental conditions. The methods were examined 
critically from the point of view their (a) experimental procedures, (b) techniques of 
parameter identification and (c) output results. 

2.3 Transient test methods for flat plate solar 
collector 

Transient methods for testing flat-plate collector can be classified into the following six 
groups based on their similarities in approach (1) simple methods, (2) indoor methods, 
(3) multi-node methods, (4) multi-test methods, (5) response function methods and 
(6) unvalidated methods 
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(a) Simple methods 

Arranovitch (1977 ) developed methods that simply added the thermal capacity term 
of the collector to its steady state energy balance. The performance equation, in terms 
of mean fluid temperature was thus given by 


q, = F'{Ta)jT - F'UiiTjm - Ta) - {rnc)^ (2.1) 

The thermal efficiency of the collector was calculated as 


qu {P^){dTjrnl dt) pif \ (^/m ~ Tfl) 

= = (ra)e - t Ui 


The efficiency was plotted against the term to give the collector param- 

eters. The method did not require controlling the inlet fluid temperature, but an 
additional experiment was necessary to determine the value of (me). The following 
limiting condition was usually considered : 


< 12°Ch-^ ( 2 . 3 ) 

at 

in order to keep the term relatively small compared to q^. The method 

was thus applicable to slowly varying insolation for which eqn (2.3) could be satisfied. 
This was clearly a limitation on the use of this method; consequently , it did not have 
universal applicability. 

Perers (1993) described a method which considered: (a) a one-node capacitance 
term to include the dynamic effects and (b) separate incidence angle modifiers for direct 
and diffuse radiation. The technique of multiple regression analysis was suggested for 
estimation of collector parameters. The methods were implemented in the testing of 
an 11 m^ long ground based (LGB) flat-plate collector and the results for collector 
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values for optical efficiencies and heat loss coefficients. This was attributed to the 
inaccurate expression for incidence angle modifiers and an incomplete second order 
model (Perers et al. ,1990; Perers and walletun, 1991). Further, scrutiny of the results 
presented by Perers and Walletun (1991) showed significant fluctuations in the values 
for effective heat capacity terms from different test runs. 

Souproun (1992) proposed a transient procedure for testing a flat-plate collector. 
It was based mainly on the analysis of a single experimental curve consisting of a 
heating and subsequent cooling process. This curve was obtained from long heating 
period of approximately 90 min or until temperature of the heat transfer fluid exceeded 
60°C, whichever was achieved first. Then the cooling process was executed by shielding 
off the sun. During both heating and cooling processes, solar radiation and ambient 
conditions remained same. Analyzing this curve, the heat losses and thermal efficiency 
of the collector were calculated as functions of the ratio of the heating to cooling times 
required to cause a temperature change between certain levels. 

A sufficiently long period of steady solar insolation was required to conduct the 
heating process. But, such steady periods occur infrequently during variable weather. 
Thus, the method was not recommended for outdoor testing. It could be adopted only 
if a solar simulator was available. 

(b) Multi-node methods 

Wijeysundera and Hawlader (1984) and Hawlader and Wijeysunder (1987) proposed 
a two-node model for the collector in which the cover capacitance was lumped with 
the plate while the fluid was treated separately. The Laplace transform technique 
was used for obtaining expressions for plate and fluid temperatures. The experiment 
was carried out by shielding the collector from the sun and suddenly circulating a 
stream of hot fluid through it. The variation of the temperature of the plate (mid 
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point) and fluid (outlet) with time was monitored continuously. With an assumed 
set of collector parameters, theoretical results were obtained and compared with the 
experimental measurements. The procedure was continued until the best agreement 
with the experimental measurements was obtained. 

With this method, no different restrictions were placed on the assumed values of 
the parameters; consequently, a large number of iterations were required. There was no 
guarantee that all guessed values were correct; different combinations of guessed values 
could lead to the same results and hence method was not dependable. Further, the au- 
thor’s results showed that a reasonable agreement was obtained for plate temperature, 
but fluid outlet temperature differs substantially from the experimental measurements. 
For such experiments high accuracy of measurement of temperature rise was needed. 

Kamminga (1985a) developed a four node model (cover, plate, fluid and insulation) 
to describe the transient behavior of flat-plate solar collectors. An approximate solution 
was obtained for the fluid temperature within the collector. In a later work (Kamminga, 
1985b), the model was simplified by considering three nodes (cover, plate k, fluid) .Finite 
Fourier transformation was applied to time varying differential equations leading to 
relations for collector characteristics. Using a series of experimental measurements 
at a fixed fluid inlet temperature, the collector parameters were determined by least 
squares methods. 

According to the author, the method led to error if Fourier variables were not chosen 
properly. Only a limited range of time interval and Fourier variables led to appropriate 
equation which were amenable to inverse Fourier transformation. Moreover results 
from different time intervals showed a relatively large scatter. 



CHAPTER 2. LITERATURE SURVEY 


15 


(c) Multi-test Methods 

Saunier and Chungpaibulpatana (1983) and saunier et al. (1985) proposed a method in 
which the solar collector was connected to an auxiliary unit and was tested under zero 
efficiency conditions. They assumed a quadratic expression for collector heat loss. The 
method required day-time test during which the collector was exposed to sun, and long 
night-time test run when several heating and cooling cycles were applied to collector. 
Both test runs were conducted outdoor under natural environmental conditions. The 
collector performance parameters were determined using multiple nonlinear regressive 
technique for experimental data. 

The method required (i) an additional test to determine the auxiliary unit parame- 
ters and (ii) a long, inconvenient night-time test for estimating effective heat capacity 
and heat loss coefficient. Moreover, according to the authors, there was significant fluc- 
tuation in the values of collector parameters obtained from different sets of test runs 
due to second order heat capacitance effect neglected in the one-node model (Saunier 
and Chungpaibulpatana ,1983). 

Chungpaibulpatana and Exell (1988, 1990) modified above method to eliminate 
the night-time test. While the theoretical basis remains the same, the experimental 
procedure was changed to determine each collector parameter individually. In addition 
to the day-time test of saunier, two additional shading tests were introduced and a 
short term heating process replaced the extended heating processes. However no ex- 
perimental verification of this method has been reported in the literature. On the other 
hand, too many tests were proposed in this method and consequently a long chain of 
data reduction step needed to be followed for estimating successive approximations of 
various terms used in calculating the collector parameters. 
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(d) Response function methods 

Rogers (1981) outlined a method based on collector’s dynamic response to actual 
weather conditions. This method enabled one to obtain the steady state character- 
istics of a flat-plate solar collector from an analysis of the experimental data. The 
method was adopted by British Standard Institution as a standard for transient test- 
ing [British Standard 6757 (British Standard, 1986)]. The models assumed that the 
thermal output of a collector during any time interval was cumulative effect of the 
fractions of energy gained during that interval and a few preceding intervals. It used 
an impulse response function to apply a proper weighting factor to solar irradiance 
incident on collector . Emery and Rogers (1984) implemented this method for toting 
three different collectors and reported that the value obtained for collector parameters 
were slightly lower than those given by steady state method. 

The method employed an experimental approach similar to that of the ASHARE OS- 
SA (ASHARE, 1986) standard and hence was much simpler compared to other methods. 
However, it sufiered from the following disadvantages: (a) it required the handling of 
large data sets (b) consequently, the amount of computation was extensive (c) the 
response function of the collector was obtained indirectly which resulted in the best fit 
to the acquired experimental data, (d) the method was not used extensively so fax. 

Wang et al. (1987) proposed a transient procedure called the filter method for 
rating the thermal performance of solar collectors. The model described the dynamic 
characteristics of the collector. The response function of the collector was determined 
from a separate shielding test. Based on this function, a recursive digital filter was 
constructed for data processing and analyzing. Although the method incorporated a 
simple experimental procedure, it required familiarity with digital signal analysis and 
filtering. 
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(f) Unvalidated method 

Fried (1990) proposed two-node and three-node collector models for transient test pro- 
cedures. The exact solutions of the models were obtained analytically to express the 
collector useful heat gain in terms of its parameters and other response functions. 
Regression analysis was suggested to estimate collector parameters and elements of 
response functions. So far, these models were not implemented in any experiments. 



Chapter 3 


PROBLEM FORMULATION AND 
SOLUTION METHODOLOGY 

3.1 Solar radiation 

Normally measuring instruments give the values of solar radiation falling on a horizontal 
surface. However, most solar equipment^ (e.g. flat-plate collectors) for absorbing solar 
radiation, are tilted at an angle to the horizontal. It therefore, becomes necessary 
to calculate the flux which falls on a tilted surface. Total amount of hourly global 
radiation reaching the tilted surface top cover is given by Sukhatme [31]. 

It = IbRb + IdRd + + Id)Rr (3-1) 

where, 

h : hourly beam radiation 
Id ■ hourly diffuse radiation 
Rb ; tilt factor for beam radiation = 

Rd : tilt factor for diffuse radiation = 
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Rr\ tilt factor for reflected radiation = 

Ib and Id are taken from the solar radiation data for India published by Mani [19]. 
Total amount of heat flux absorbed by the absorber plate is given by 

S = hRb {Ta)b + ( {hRd) + {h + Id) Rr) {ra)d (3.2) 

where, 

{Ta)b : transmissivity-absorptivity product for beam radiation 
(rQ:)d : transmissivity- absorptivity product for diffuse radiation 

3.2 Governing differential equations and boundary 
conditions 

This section provides a more appropriate and accurate multidimensional transient anal- 
ysis for double glazed flat plate solar collector (see Fig. 1.1). Following assumptions 
are made in the analysis : 

1. Temperature drop through the glass covers is negligible. 

2. The sky is considered as a black body. 

3. Variation in all properties with respect to temperature is neglected. 

4. Accumulation of dust and dirt on the collector is disregarded. 

5. Shading of collector plate is not considered. 

6. Temperature on periphery of the pipe is considered to be same as that of the 
absorber plate. 
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3.2.1 Cover2 

The cover thickness is so small that it is reasonable to consider a uniform temperature 
throughout. (The lumped heat capacity approach also yields, Bi < 0.1). Hence, we 
get a two dimensional distribution of temperature , Tc 2 {x, y), in cover2. The governing 
equation can be obtained from an energy balance in a small volume of thickness, Sc 2 , 
of the cover2 as, 


fcc2l“?r^ + 


dx^ 


dy^ 


, Qa-c2 , QcI—c 2 , ^ dTc2 

+ ! h Qabs = Pc2Cc2' 


^c2 


5c2 


dt 


(3.3) 


where, 

Pc 2 : density of cover2 

Cc 2 • specific heat of cover2 

kc 2 : thermal conductivity of cover2 

qa-c 2 '■ heat flux from ambient to cover2 

= convective heat flux ( from ambient to cover2 ) + 

radiative heat flux ( from ambient to cover2 ) 


— hwiTa — Tc2m) + <^^c2{T^ky ~ ^c2m) 
qci-c 2 ■ heat flux from coverl to cover2 
= average convective heat flux ( from coverl to cover2 ) 4- 
average radiative heat flux ( from coverl to cover2 ) 

= hci-c2{Tclm — Tc2m) + 

‘c2 ' 

[l—ezpC— — — ^)] 

qabs • heat absorbed by cover2 per unit volume = L - 

The symbols Ta, Tdm, Tc 2 m, ^ci, ^c 2 and K stand for ambient temperature, mean 
coverl temperature, mean cover2 temperature, emissivity of coverl, emissivity of cover2 
and extinction coefficient of cover glass 2 respectively. An empirical expression for con- 
vective heat transfer coefficient between cover2 and ambient, h^,, proposed by Sparrow 
and Tien has been used [29j. The symbol hd_c 2 stands for convective heat transfer 
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coefficient between coverl and cover2 and is obtained by using the empirical expression 
suggested by Buchberg et al. [4]. The heat absorbed by cover, qai*, is calculated using 
Bouger’s law for a partially transparent medium |9]. 

It may be noted that the thermal radiation term from ambient to cover is calculated 
considering the ambient air as a black body which emits at sky temperature. The sky 
temperature is calculated from empirical relation given below [31] : 

Tsky = Ta - 6 

Boundary condition is taken as Adiabatic (Neumann type), i.e. : 

^ = 0 (3.4) 


3.2.2 Coverl 

As in the case of the cover2, here also we consider a two dimensional variation in the 
temperature of the coverl i.e., Tci{x,y). The governing equation can be obtained from 
an energy balance in a small volume of thickness, dci, of the coverl as, * 


kcil-^TTr + 




9y2 


, 9c2-ci , 9p-ci , dTci 

T J 1 J r Qabs — 


-'cl 


^cl 


where, 

pci ■ density of coverl 

Cci : specific heat of coverl 

kci : thermal conductivity of coverl 

qc 2 -ci ■ heat flux from cover2 to coverl 

= average convective heat flux ( from cover2 to coverl ) + 

average radiative heat flux ( from cover2 to coverl ) 

= hcl-c2{Tc2m - Tclm) + 

qp-c 2 ■ heat flux from plate to the coverl 


(3.5) 
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= convective heat flux ( from plate to coverl ) + 

radiative heat flux ( from plate to coverl ) 

L ’’P _L 

— J-clm) * TIlIIjrXrZTT" 

1 1 — **** ^1 

Qabs ■ heat absorbed by coverl per unit volume = L 

The symbols Tdmj Tc 2 m, Tpm, €ci> and K stand for mean coverl temperature, 
mean cover2 temperature, mean temperature for absorber plate ,emissivity of coverl, 
emissivity of plate and extinction coeSicient of cover glass 1 respectively. The symbol 
hp_ci stands for convective heat transfer coefficient between plate and coverl and is 
obtained by using the empirical expression suggested by Buchberg et al. [4]. 

Boundary condition is taken as adiabatic (Neumann type), i.e. : 


dTci 

dn 


= 0 


(3.6) 


3.2.3 Absorber plate 

.A.S in the case of the covers, here also we consider a two dimensional variation in 
the temperature of the absorber plate i.e., Tp{x^y). The governing equation can be 
obtained from an energy balance in a small volume of thickness, 6p, of the plate as. 


kl^ 


+ 


5% 

5y2 


+ 


Qci-p 

Sp 


+%2+ 


if-p I ^ 


where, 

Pp : density of the absorber plate 

Cp : specific heat of the absorber plate 

kp : thermal conductivity of the absorber plate 

S : solar flux absorbed by the absorber plate 

Qp-ci '■ heat flux from coverl to the plate 

= average convective heat flux ( from coverl to plate ) + 


(3.7) 
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average radiative heat flux ( from coverl to plate ) 

__ u (rp rp \ . cr(T^j^— 

— "p-cU-tclm ~ J-pm) + f J_rx_i^ 

^‘ci ‘P ^ 

Qi^p : heat flux from insulation to absorber plate = —ki^ |z=ii 
qf-p : heat flux from fluid to absorber plate = hf{Tf ~ Tp) 

Heat transfer coefficient between plate and fluid, h /, is calculated using the assump- 
tion of fully developed laminar flow through the pipe and constant heat flux condition 
|15|. 

Boundary condition is taken as Adiabatic (Neumann type), i.e : 


3.2.4 Insulation 

The temperature distribution Ti{x,y,z) in the insulation is considered to be three 
dimensional ( Lumped heat capacity method cannot be applied as, B,- > 0.1 ). The 
governing differential equation is obtained from an energy balance in a small volume 
of thickness Si, of the insulation as. 


, ^d^Ti , d% . d^Ti^_ ^^dTi 

5x2 + Qy2 + 5^2 1 Qi 


(3.9) 


where. 

Pi : density of the insulation 

Ci : specific heat of the insulation 

ki : thermal conductivity of the insulation 

There is no heat generation or absorption inside the insulation. The boundary 
conditions are; 


Faces exposed to the atmosphere (Robin’s type B.C.) : 
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dT- 

-H-^) = hUT,-Ta) 


(3.10) 


Face in contact with the absorber plate (Dirichlet type B.C.) 


T = T 

J-l — J.p 


(3.11) 


3.2.5 Fluid 

The fluid flow in the duct is considered to be one dimensional i.e. average values of 
temperature, T/, are used in each cross sectional area of the fluid stream. The governing 
differential equation of the fluid flow is : 


where, 

Pf : density of the fluid 
Cf : specific heat of the fluid 
m : mass flow rate of the fluid 

Qp-f : convective heat flux from the absorber plate to the fluid 

= />;(r,-r,) 

P = perimeter — itDi 
Af = area of cross-section = 

The inlet boundary condition for the fluid at x=0 , y = 0 is : 


(3.12) 


T, = Tfl 


(3.13) 
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3.3 Method of solution 

The finite difference method is used to discretize the governing differential equations. 
The technique used for discretization is Alternative Direction Implicit (ADI) for cover2, 
coverl, absorber plate, insulation & Implicit scheme for fluid. In the method, the entire 
region is divided into small grids of finite size. Temperatures are calculated on all the 
nodal points at each time step as we move forward in time domain. In the nomenclature 
used, i, j and k represent parameters in x,y and z directions respectively and n is the 
parameter in time domain. The discretized elements of covers, plate, insulation and 
fluid of a solar collector are shown in Fig. 3.1. The computational algorithm for 
the coupling solution is indicated by a flowchart (Fig. 3.2). Partial derivatives are 
expressed in the following forms : 


dx 


Tj+i — Tj-i 

2Ai 


+ 0{Ax'^) 


dx"^ 


71+i-2r,+T,-i 

2Ax2 


+ 0(Ar2) 


dt 




At 


+ 0(At) 


The notation 0(Aa:^) means that arriving at the above, order of (Ax^) and higher 
have been neglected. 0(Aa;^) is called truncation error. Truncation error for 2-D 
ADI is 0(At^, Az^, Ay^) and for 3-D ADI it is (At,Ax^,Ap^,Az^) [32j. Truncation 
error is the difference between the exact mathematical expression and the numerical 


expression. 
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Ta,j+i) 

✓ 



(a) Discrete grid points in covers and plate 


T;(i,j,k+1) 

<?> 



(b) Discrete grid points in the insulation 


T(i.l) 

e 


(c) Discrete grid points in fluid 



T(i+1) 

O 


T(i+l,j) 


. T(i,j+l,k) 

e---T(i+l,j,k) 


Figure 3.1: Discrete Grid Points of Different Components of a Solar Collector 
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3.3.1 Cover2 : 

A set of discretized equations is obtained by numerical approximation of the governing 
equation (3.3) and its boundary condition (3.4) into each one of the control volumes 
of the grid. The evaluation of the heat conducted is made using Alternative Direc- 
tion Implicit (ADI) technique. The convection heat exchanged with outside is taken 
implicitly and remaining heat exchanges are taken explicitly. 

Two equations are established for each point of the grid, each one implicit in one 
direction. 

1. First implicit equation in x direction (move ^ forward in time domain) yields 


bTc 2 {i - l,i,n + 1/2) + dTcaihj, n -t- 1/2) + aTc 2 {i + 1, j,n + 1/2) = r (3.14) 

where, 
b = kc2^y^ 

d = -[2kc2^y^ + 2pc2C'c2AxAy|| -f hyjAxAy] 

a = ^c 2 Ay|| 

r = -kc 2 Ax^[T^{i,j -l,n) - 2 Tc 2 (i,j, n) + ^c 2 (^,i 4- 1, n) - AxAyl/i^uTa 
-l-aec2(r/fcj, - Tjs(z,j,n))] - AxAy[hci.c2(Tcim - T^m) + 

'‘«cl «c2 ' 

-qahsAxAy5c2 - pdiCc2AxAy5c2 

2. Second implicit equation in y direction (move by another— in time domain) 
yields 


bTaihj - 1, n -}- 1) -1- dTc 2 {i,j,n + 1) 4- aTc 2 {iJ 4- 1, n 4- 1) = r 
where, 

b = kc2Ax^ 


(3.15) 
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d — — ‘^Pc2Gc2^^Aky + /lyjAxAy] 
a = 

r = — /Cc 2 Aj/^|[Tc 2 (i — + 1/2) — 2Ta{i,j,n + 1/2) + Tc 2 {i,j + l,n + 1/2) 

-AxAy[hyjTa + - T^2ihj,n+ 1 / 2 ))] - AxAy[hci-a{Tcim - Tc2m)+ 

~ 9a6sAa:Ay5c2 - Pc2Cc2 Ax Ay(5c2 

The coefficient matrices of the resultant equations are tridiagonal. The coefficients 
a, d, b are respectively for upper diagonal, diagonal, lower diagonal elements and r 
is a constant vector. Therefore, it can be solved by a Tridiagonal Matrix Algorithm 
(TDMA). The boundary conditions in discretized form can be expressed as (by image 
point boundary): at 


i = 1 

Tc2{i — 1) j) = Tc2{i + 1, j) 

i = rix 

Tc2{i + 1) j) = Tc2{i — l,i) 

3 = 1 

Tc2{hj — 1) = Tc2{i,j + 1) 

j = ny 

Tc2ii,j + 1) = Tc2{i,j - 1) 

3.3.2 

Coverl 


A set of discretized equations is obtained by numerical approximation of the governing 
equation (3.5) and its boundary condition (3.6) into each one of the control volumes 
of the grid. The evaluation of the heat conducted is made using ADI technique. Heat 
exchanges on cover2 side and the plate side are calculated explicitly. 

Two equations are established for each point of the grid, each one implicit in one 
direction. 

1. First implicit equation in x direction (move ^ forward in time domain) yields 
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bTci{i - l,j,n+ 1/2) + dTci{i,j,n + 1/2) +aTci{i+l,j,n+ 1/2) = r (3.16) 


where, 



d = —[2kciAy^ + 2pciC'ciAxAj^^ 



r — j 1)?^’) 27).i (i, j, Ti) -|- Tgi (z, j 1, n) AxAj/[/ici— c 2 

(Tc2m - Tcim) + “ AxAy[hp-cl{Tpm “ rci,„) + 

'‘«cl «c2 

‘"/feiS' r “ 9a6.AxAy(5ci - 2pciC'c, AxAt/(5ci^^i^ 

^<cl «P ' 

2. Second implicit equation in y direction (move by another^ in time domain) 


yields 


bTci{i,j — 1, n + 1) + dTci{i, j,n + 1) + aTci{i,j + 1, n + 1) = r (3 A?) 

where, 

b = kciAx^ 

d = — [2fcciAx^ + 2pciC'ciAxA7/^] 

a — kciAx-^ 

r = -kciAy^[Tci{i - l,j,n + l/2) - 2Tci{i,j,n + 1/2) + TdiiJ + l,n+ 1/2)] 

-AxAy[/ici_c2(3;2„» - Tcim) + - AxAy[hp_ci(rp„ - Tam) 

'^«cl «c2 ^ 

- <i-bsAxAy5a - 2paCaAxAy5a'^^^^^^^^ 

The coeflScient matrices of the resultant equations axe tridiagonal. The coefficients 


a, d, b are respectively for upper diagonal, diagonal, lower diagonal elements and r 
is a constant vector. Therefore, it can be solved by a Tridiagonal Matrix Algorithm 
(TDM A). The boundary conditions in discretized form can be expressed as (by image 
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point boundary): at 
i = l T,i(i-l,j)=T,,{i + l,j) 

i = nx Tci{i + lJ) =Tci{i- l,j) 

j = 1 Tci{i,j - 1) = Tcx{i,j + 1) 

j = Uy Tci {i,j + 1) = Tci {i, j - 1) 

3.3.3 Absorber plate 

The governing energy equation for absorber plate is (3.7). A discretized equation is 
obtained by numerical approximation of the governing equation, using ADI criteria. 
Heat exchanges on cover side and with fluid and insulation are calculated explicitly. 

Two equations are established for each point of the grid, each one implicit in one 
direction. 

1 . First implicit equation in x direction (move ^ forward in time domain) yields 


bTp{i- I, j,n + 1/2) + dTp{i,j,n + 1/2) + aTp{i + l,j,n + 1/2) = r 


(3.18) 


where, 
b = kpAy-^ 

d = -[2kpAy^ + 2ppCpAxAy 
a = kpAy^ 


Afj 


= -/cpAa;^[Tp(i, j - l,n) - 2Tp{i,j,n) + Tp{i,j + l,n) - AxAy[hp-ci{Tcim 

■Tp^) + "I'i.'lLhh - - l.n) 

-hfPAy[Tf{i, 3 ,n) - Tp{i, 3 ,n)] - SAxAy - 2ppCpAxAy6p^^^ 
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2. Second implicit equation in y direction (move by another^ in time domain) 
yields 


bTp(i,j - l,n + 1) + dTp{i,j,n + 1) + aTp{i,j 4- l,n+ 1) = r (3.19) 

where, 
b = kpAx-^ 

d = —[2kpAx-^ + 2ppCpAxAy^] 

r = —kpAy^[Tp{i — l,j,n + 1/2) — 2Tp{i,j,n + 1/2) + Tp{i,j + l,n + 1/2)] 
-AxAy[/ip_ci(Tcim - Tpm) + - ki^^[Ti{i,j,n 2 - l,n) 

-Ti{i,j, n^,n)]-hfPAy[Tf{i,j,n)-Tp{i,j,n)]-SAxAy-2ppCpAxAy6p'^^^^^^^^f^ 
The coefficient matrices of the resultant equations are tridiagonal. The coefficients 
a, d, b are respectively for upper diagonal, diagonal, lower diagonal elements and r is 
a constant vector. It can be solved by a Tridiagonal Matrix Algorithm (TDMA). The 
boundary conditions in discretized form can be expressed as (by image point boundary): 
at 

i = 1 7],(z — 1, j) = Tp[i + !,_)) 

i = Tp{i + 1, j) = Tp{i - l,j) 

j = l Tp{i,j-l)=Tp{i,j + l) 

j = Uy Tp{i,j + 1) = Tp{i,j — 1) 

3.3.4 Insulation 

The discretization is carried out by means of a regular three dimensional grid. For each 
one of the control volume elements, a discretized equation is obtained as a result of 
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numerical approximation of the governing equation (3.9) and its boundary conditions 
(3.10 & 3.11). ADI method is used for evaluation of the heat conduction. Three 
equations are obtained for each point of the grid, each one implicit in one direction. 

1. First implicit equation in x direction (move ^forward in time domain) yields 


bTi{i — l,j,k,n + 1/3) + dTi{i,j,k,n+ 1/3) + aTi{i + l,j,k,n + 1/3) = r (3.20) 
where, 

0 — Ax 

d = -[2ki^ + 3piCi^^] 

° Ax 

r = - l,k,n) - 2Ti(iJ,k,n)+T,{i,j + l.fc.n)] - 

{Ti{i,j,k-l,n) - 2ri(i,j,A;,n) -Ti(i,j,(: + l,n)] 

2. Second implicit equation in y direction (move by another^ in time domain) 
yields 


bTi{i,j — l,/c,n + 2/3) + dTi{i,j, k,n + 2/3) +aTi{i,j + l,fc,n + 2/3) = r (3.21) 


where, 

L _ jl.AiAz 


a = ki 


AiAz 

As/ 


r = -ki^^[Ti{i — lj,k,n+ 1/3) - 2Ti{i,j,k,n + 1/3) + Ti{i + l,j,k,n+ 1/3)] 
-ki^[Ti{i,j, k - l,n+ 1/3) - 2Tiii,jAn + 1/3) - Ti{i,j,k + l,n + 1/3)] 
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3.Third implicit equation in z direction (move by another^ in time domain) yields 


bTi{i,j,k- l,n + 1) + dTi{i,j,k,n+ 1) + aTi{i,j, k + l,n+l) = r (3.22) 

where, 

d=-l2k,^+Zp,C,^^] 

_ L . AxAy 
— >^1 Az 

^ = -ki^[Ti{i - l,j, k,n + 2/3) - 2Ti{i,j, k,n + 2/3) + Ti{i + l,j, k,n + 2/3)] 
j-l,k,n + 2/3) - 2Ti{i, j, k,n + 2/3) - Ti{i, j + l,k,n + 2/3)] 

The coefficient matrices of the resultant equations are tridiagonaJ. The coefficient 
a, d, b are respectively for upper diagonal, diagonal, lower diagonal elements and r is 
a constant vector. It can be solved one by one using Tridiagonal Matrix Algorithm 
(TDM A). The boundary conditions are : at 

i = 1 Ti(i - l,j,k) = Ti(i + l,j.k) + -T.l 

i = n, Ti[i + l,i,k) = T,(i -l,j,k)- ^[T,{i,j, k) - rj 

j = 1 T,(i,j - 1,A:) = T|(i,i + l,k) + -T.l 

j = + l,fe) = Ti(i,i -l,k)- 5^[T,(i, j,i) -T.l 

k = 1 k-l)= T,(iJ, k+l) + k) - T.) 

k = n^ Ti(i,j,k) =Tp(i,j) 
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3.3.5 Fluid 

For each one of the finite control volumes that results by subdividing the pipe, the 
discretization equations are obtained by numerical approximation of energy equation 
(3.12). The internal energy transportation is calculated implicitly. 

lmC, + p,A,C,j^ + h,PAx] 

(3.23) 

The boundary condition is ; at i=0, j=0 

^/(*1 j ) ~ 

Equations described above yield temperature distributions in cover2, cover 1, ab- 
sorber plate, insulation and fiuid as a function of time and spatial co-ordinate. The 
temperature distribution can be found out at time t-l-At if we know the temperature 
distribution at time t. The initial temperature of different elements of the collector 
is taken as the ambient temperature when the calculation is initialized. The accuracy 
improves when the spatial and temporal discretization densities are increased. The nu- 
merical solution tends towards an asymptotic value when finer grids are used. For the 
discretization used, further increment in the density of discretization does not cause 
appreciable improvement in the numerical solution. Two important parameters namely 
the stagnation temperature and efficiency' of the collector are defined below'. 

3.4 Stagnation temperature 

It refers to the condition of a fiat plate solar collector when the liquid fiow through the 
collector is stopped. There is no useful heat gain and its efficiency is zero. In this case. 
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absorber plate attains a temperature such that the energy absorbed by it is equal to 
the heat loss. This temperature is the highest that the absorber plate can attain and 
is called the stagnation temperature. 

3.5 Efficiency 

It is an important parameter to estimate the performance of a collector. The instanta- 
neous efficiency of a flat plate solar collector is defined as the ratio of useful heat gain 
to the flux incident on the collector. The mean efficiency of collector is defined as the 
ratio of the total useful heat gain to the total incident radiation on the collector over 
the period of its operation. The collector efficiency based on the absorber plate area is 

_ ThCpjTfo - Tfi) 

^ Apit 
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Figure 3.2: Computational Algorithm of Numerical Solution 















Chapter 4 


EXPERIMENTAL SET-UP AND 
TEST PROCEDURE 

4.1 Purpose of experiment 

The experiments were conducted with the purpose of getting real time data for valida- 
tion of the theoretical model developed here for the flat-plate solar collector. 

4.2 Description of experimental set-up 

The experimental set-up was completed with the following major components. 

1. Absorber plate 

2. Cover plate 

3. Insulation 

4. Casing 

5. Instrumentation 


37 
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1. Absorber plate 

The choice of material depends on the thermal properties (conductivity & capacitance), 
mechanical properties (density, ductility & ease of soldering/welding), corrosion resis- 
tance, cost and availability . Galvanized iron, steel, aluminum, copper, stainless steel 
and plastic etc. may be used as material for plate. In the present study, copper was 
used. Although it is expensive, but it has good thermal and mechanical properties and 
also corrosion resistance. 

Three different types of absorber plate configurations normally used are (1) Pipe 
and fin (2) Full water sandwich (3) Semi water sandwich type as shown in Fig. 4.1. Pipe 
and fin type (all copper) configuration was taken due to its long life, low maintenance 
and capital cost, good corrosion resistance and its ability to withstand high pressure . 
The absorber used in present test rig is shown in Fig. 4.2. 


— o — o — o — o~ 

(a) pipe and fin type absorber plate 



(b) Full sandwich type absorber plate 


(c) Semi - sandwich type absorber plate 


Figure 4.1: Absorber Plate Configuration 
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2. Cover plate 

Glass is the most commonly used cover material. It has high transmittance (of the 
order of 0.85 or more) in short wavelength range and is practically opaque to long wave 
infrared radiation emitted by the absorber plate. Glass deterioration is negligible even 
over very long periods of exposure to intense ultraviolet radiation. The cover plate 
used for collector was made of simple window glass of 4 mm thickness. 

3. Insulation 

Several types of insulation materials are used in collectors to minimize the heat loss 
from the back and sides. The commonly used materials are fiberglass, ceramic fiber 
blanket, mineral fiber blanket, calcium silicate foam, urea formaldehyde foam, poly- 
urethane foam and glass wool. Glass wool was chosen in the present case since it was 
inexpensive and easily available. 

4. Casing 

The material chosen for casing was a GI sheet of 30 gauge thickness because of good 
strength, good corrosion resistance, low cost and easy availability. 

A complete assembly of the collector is shown in Fig. 4.3 and a line diagram of the 
complete test rig is shown in Fig. 4.4. 

5. Instrumentation 

(a) Temperature measurement 

The temperature measurement was done using 24 gauge copper constantan thermo- 
couples (supplied by Love Controls Corpn. USA ). The location of thermocouple is 
shown in Fig. 4.2. The temperatures were recorded using direct temperature measur- 
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mg instrument. The thermocouples were calibrated using constant temperature water 
bath between 0-C and 100"C. Thermometer waa used to meaaure the ambient, the 

inlet and the outlet water temperatures. The accuracy of temperature measurement 
was ±0.2°C. 
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Water tank 



(b) Mass flow rate measurement 

Mass flow rate was measured by direct measurement of volume discharged, with the 
help of a calibrated measuring flask. The time was taken by a stop watch (made by 
Hanhart, Germany ). The accuracy of volume measurement was ±5 ml and that of 
time was ±0.1 s. 


4.3 Fabrication and installation 

Generally the shape of a solar flat-plate collector is rectangular. The fabrication of 
the collector was done in several stages. The main structure was started by making 
a rectangular box of wood of overall dimension 1.05 mxl.05 m. Further the absorber 
plate was fixed in the structure by giving space for insulation. The insulation used in 
flat-plate collector was glass wool. The thickness of insulation used was 80 mm for the 
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back and 40 mm for the sides. 

The absorber plate was made of copper material of size 1.05 mx 1.0 m. To get 
better surface contact between plate and pipe, nine depressions made in the plate at a 
pitch of 85 mm, Nine copper pipes were placed in the depression and were brazed with 
the plate filling the depression with copper eutectic. The copper pipes were taken of 
12.5 mm (1/2 inch) in diameter and these were brazed in series using elbows. They 
were then tested for leakage using the tap water. The upper surface of the absorber 
plate was painted with black board black paint. Five Thermocouples were provided at 
four corners and one in the middle. The final size of collector plate so obtained was 
1.0 mx.95 m. Provision was made for the movement of absorber plate by screw and 
pinion arrangement for varying the distance between the absorber plate and the cover. 

The cover plate was provided above the absorber plate which helped to prevent 
convective and re-radiative heat losses to the atmosphere. The option w'as available for 
using either single or double glazing. There was also provision for varying the distance 
between the cover plates in a double cover design. Finally the sides of the collector 
were covered with GI sheet of 30 gauge. 

The stand for collector was fabricated, the frame of which was made of slotted angle 
iron. The bearing type arrangement shown in fig. 4.3 was made for the movement of 
the flat plate collector. 

Once the solar flat-plate collector was ready. A PVC tank of 500 liter capacity was 
connected to the flat plate collector through proper fittings. A float valve was used to 
maintain a constant water level. After careful and thorough investigation, a suitable 
site of the set-up was decided. There was availability of electricity and water tap on the 
roof of heat transfer lab which was chosen for the installation. Finally, Initial testing 
was done to get the correct orientation and full sun shine. 
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4.4 Test procedure 

The experiments were conducted during the winter, months namely Dec, Jan and 
Feb on the set-up. The opening of control valve used in the pipe line for the control of 
mass flow rate through the flat-plate solar collector was regulated and the temperatures 
readings were taken at regular intervals of one hour each from 9 A.M. to 4 P.M., on 
daily basis. 

(a) Collector inlet 

(b) Collector outlet 

(c) Plate temperature at 5 different location 

(d) Ambient temperature 

The direct measurement of volume of water was taken at outlet of collector. The 
time was taken by a stopwatch for collection of water. The data were repeated for 
three different flow rates namely 15.5 kg/h, 31.5 kg/h and 56.5 kg/h. 

The experiment was conducted only on clear days to overcome the influence of 
cloudiness on the incidence of solar radiation. 


4.5 Experimental uncertainties 

Errors in the experiment may arise due to various reasons. These can be classified as: 

(a) Systematic or determinate errors, arising as a result of poor calibration of in- 
struments and 

(b) Random or unavoidable errors, arising as a result of uncertainty associated 
with the accuracy of measuring instrument and the ability of the operator to read such 
instruments both accurately and consistently. 

While no attempt was done in the present study to quantify these errors, all feasible 
efforts within the practical limits were implemented to ensure that the possible errors 
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are minimized. 


Uncertainty Analysis 


The data reduction equation of efficiency is given by 

_ mCp{Tio-Tfi) 

” “ Aph 

The value of m, T/o and T/j were measured while the values of Ap and It were taken 
directly. Therefore, the efficiency depended only on rh, T/o and T/j. 
i.e. T) = r){Tfi^ Tf^, m) 

The uncertainty in T/t and T/o are 0.2°C. The uncertainty in rh can be determined 
as follows : 


7n = m(V,t) 

Where both V and t were measured. The uncertainty in V and t were 5 ml and 0.1 
s. The uncertainty analysis expression in the measurement of m can be written as [7] : 

= (^F + (^F 

— (';005'i2 4_ 

“ v .25 ''' V58/ 

= .000403 


(^) = .020007 

Uncertainty in the measurement of mass flow rate was 2.0%. The uncertainty 
analysis expression in the efficiency can be written as [7] : 

= {W+iW+im'^ 

= .0005814 


(^) = .02411 

The uncertainty in efficiency was 2.41% at 12 noon for mass flow rate 15.5 kg/h. 
On a particular day for a given mass flow rate, the uncertainty will vary from 
morning to evening because of changing solar insolation and therefore the temperature 
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rise. Hence one can talk about some minimum or maximum value in this situation. 
This uncertainty falls in the range of 2.41% to 7.35%. 

At a given time of day for different mass flow rate, the uncertainty will be different 
on temperature rise. This value falls in the range of 2.41% to 4.51% at noon for mass 
flow rate in the range of 15.5 kg/h to 56.5 kg/h respectively. 



Chapter 5 


RESULTS AND DISCUSSION 

5.1 Introduction 

The governing equations for transient modeling of a fiat plate solar collector were 
developed in section 3.2. These were discretized and the complete set of discretized 
equations along with the B.C.’s (3.14 ) through (3.23) were solved by finite difference 
technique. The computed results are presented in the following sections and are com- 
pared with the experimental data collected during the transient testing of the solar 
collector. 

5.2 Results 

The performance of fiat plate solar collector for a given day, location and position 
depends on inputs namely mass fiow rate of fluid, inlet fiuid temperature, starting 
time, closing time and initial temperature distribution in different components of the 
collector. 

The experiments in the present work were conducted during the winter, months 
namely Dec, Jan, and Feb of an year at Kanpur (India) as discussed earlier. The 


46 



CHAPTER 5. RESULTS AND DISCUSSION 


47 



9 10 11 12 13 14 15 16 17 

Time (h) 

Figure 5.1: Variation of Solar Flux on Top Cover of a Tilted Surface 

experiments were conducted from 9 A.M. to 4 P.M. The data obtained from experiments 
were compared with the computed results. In the experiment, the inlet temperature of 
the fluid and ambient temperature were continuously varying which also was considered 
in the theoretical model. For computations, the initial temperature distribution of 
different components of a flat plate collector was assumed to be same i.e. equal to 
the ambient temperature. The dimensions of collector, material properties and other 
parameter values used in computations are given in ’Appendix A’. 

The variation of incident flux on the top cover of tilted flat plate collector is shown 
in Figure 5.1. Figures 5.2 and 5.3 show the comparison of the theoretical model results 
with the experimental results for single cover and double cover respectively. It gives 
the variation of the rise in the fluid temperature with the time of the day. For the 
single cover the maximum temperature rise is about 20°C and for double cover it is 
close to 22°C. The maximum temperature rise occurs between 12 noon and 1 P.M. 
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when the solar radiation flux is maximum. The above figures show good agreement 
between the experimental and the theoretical results from 9 A.M. to 2 P.M., but there 
is discrepancy during 3 P.M. to 4 P.M. This could be attributed to various factors a 
few of which are given below ; 

1 . The solar insolation taken in the theoretical model may be higher than the actual 
values. 

2. The effect of dust is not considered in the theoretical model which may affect 
more significantly towards the end of the day. 

Figures 5.4 and 5.5 show the variation of plate mean temperature with the time of 
the day for single and double covers. The maximum plate temperatures are about 318 
K and 319 K for single and double covers respectively. The model predicted maximum 
plate mean temperature was obtained between 1 P.M. and 2 P.M. Overall the model 
predication are in good agreement with the experimental data. 

Figures 5.6 and 5.7 show the variation of instantaneous efficiency with the time of 
the day for single and double covers. There is a continuous decrease in efficiency with 
the time of the day. It is due to the increase in inlet fluid temperature while taking 
the experimental data as the day progresses and this aspect has been considered in the 
theoretical model as well. Figure 5.8 displays the nature of instantaneous efficiency 
with the time of the day when inlet temperature of fluid remained constant. It is 
seen that efficiency first rises from 9 A.M. to 11 A.M. beyond which it remains nearly 
constant till about 3 P.M. and then starts decreasing. This can be explained as follows. 
The collector is in transient state during start and much of the incident flux is absorbed 
by covers, absorber plate and insulation in raising their temperatures before they reach 
a quasi-steady state. Thus the useful energy recovered is low during this period and 
hence efficiency. So the period between 11 A.M. to 3 P.M. can be taken as the period 
during which quasi-steady state can be assumed. 
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Instantaneous Efficiency 
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Figure 5.7: Instantaneous Efficiency with Double Cover 







Instantaneous Efficiency 
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X • ^mnnrtant parameter which strongly influences 
The inlet fluid temperature is an importa p 

tke performance of a flat plate solar collector. The effect is iUastrated by plotting ^ 
„ean efficiency with the inlet fluid temperature which was varied from 20"C to 60 
(see Fig. 6.9). It is observed that efficiency drops sharply at an —g rate ( rom 
0.66 to 0.24) with increasing inlet fluid temperature (from 20-C to 60=C) as reported 

in literature [31]. . „„„„ 

, f thP experimentation used a continuous copper 
The flat plate collector fabricated for the p 

t r . increased continuously as it passed through. 

pipe. The temperature of fluid therelor 

u ni«tP is shown as a 3-D distribution in Fig. 

This temperature distribution of absor er p c+u fiovfrv 

• f ct.mation temperature with time of the day for 
5.10. Figure 5.11 shows the variation of stagnation x p ^ ^ 

• f 1 P.nlts show that stagnation temperature is maximum 
double cover. The experimental result 

U . ,hB experimental and theoretical values .s obvrons 
at 2 P.M. Good agreement between the exp 


from the figure. 
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Inlet Fluid Temperature ( C) 

Figure 5.9: Mean Efficiency with Inlet Fluid Temperature with Double Cover 


Theoretical 


Temperature (K) 



Figure 5.10: 3-D Temperature Distribution in Absorber Plate 
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Figure 5.11: Stagnation Temperature with Double Cover 


Figure 5.12 displays the rise in fluid temperature with time of the day for three 
different flow rates namely 15.5 kg/h, 31.5 kg/h and 56.5 kg/h. First the temperature 
rises till the noon, when the solar radiation flux is maximum. The maximum tempera- 
ture rises for the lower flow rates. Figure 5.13 demonstrates the variation of plate mean 
temperature with the time of the day for three different flow rates. The figures 5.12 
and 5.13 show good agreement between the experimental and theoretical values until 2 
P.M. only after which the deviation between the two increases. Figure 5.14 shows the 
variation of instantaneous efficiency for the three different flow rates. The maximum 
efficiency is obtained for the larger flow rate. Fig. 5.15 displays the variation of mean 
efficiency with the mass flow rates. From the graph it is clear that mean efficiency 
increases as the mass flow rate increases. But for flow rates greater than 40 kg/h the 
mean efficiency becomes approximately constant. 

Figures 5.16 and 5.17 show the variation of rise in temperature and plate mean 
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Figure 5.14: Instantaneous Efficiency for Different Mass Flow Rates 
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Figure 5.15: Mean Efficiency with Mass Flow rate with Double Cover 
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temperature tor four differeut tilt angles namely 12”, 24”, 36" and 48". tilt angle 
of 36" gave the maximum temperature rise and plate mean temperature for most of 
the time during the day. This suggests that a tilt angle of 36" i.e. (4> + 10) may be 
considered as an optimum as already established in the literature 131]. 

Figure 5.18 shows the variation of mean efficiency with varying the distance between 
the plate and the cover. It is obvious from the figure that as the spacing between the 
cover and the plate increases the mean efficiency also increases. Fig. 5.19 shows the 
mdation of the mean efficiency with varying the distance between the covers m a 
double cover design. The graph shows that greater the spacing between the covers 
Mgher in the mean efficiency. This increase in mean efficiency in Figures o.l8 an • 
he attributed to the decrease in the value of convective heat transfer coefficient 
,vi.h increase in spacing between the plate and cover or between the two covers m a 

double cover design [4]. 


can 

wii 
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Figure 5.17; Plate Mean Temperatiu-e for Different Tilt Angles 



Figure 5.18: Mean Efficiency with Spacing between Plate and Cover 
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Figure 5.20 shows the variation of uncertainty in efficiency with the time of day. 
Figure shows that uncertainty is minimum at 12 noon due to maximum temperature rise 
that occurs at this time. The uncertainty increases after 12 noon due to progressively 
decreasing rise in temperature. 


5.3 Conclusions 

• The rise in fluid temperature is maximum between 12 noon to 1 P.M. (close to 
12 noon). 

• The plate mean temperature is maximum between 1 P.M. to 2 P.M. (close to 2 
P.M.) 

• The mean efficiency decrease as inlet temperature of fluid increases 

• The mean efficiency is constant approximately for mass flow rates greater than 
of 40 kg/h which indicates the minimum flow rate of the collector. 

• The tilt angle 0 + 10 is optimum tilt angle as established in literature. 

• The period between 11 A.M. to 3 P.M. can be taken as the period during which 
quasi-steady state may be assumed for the solar collector. 

• The performance of the solar collector with double cover is better than that with 
single cover. 

• The spacing between the absorber plate and cover and between the two covers 
should be large. The recommended values are 50 mm to 80 mm and 50 mm 
respectively. 

• The theoretical results are in good agreement with the experimental data from 
9 A.M. to 2 P.M. which imply that the transient model developed for the flat 
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plate solar collector has been validated by the experimental results of the present 
work. This model therefore can be used to characterize a flat plate solar collator 
as a function of all the relevant parameters. 

5.4 Suggestions for future work 

• A similar transient analysis can be extended (both experimentally and theoret- 
ical) for some novel design like honeycomb collector, evacuated tube collector, 
plastic collectors and also with selective surfaces. 

• The effect of variation of physical properties with temperature of the materials 
involved. 

• The effect of shadow and dust on top of the collector can be investigated. 

• Present model considers the temperature on the periphery of the pipe same as 
that of the absorber plate i.e. the thermal resistance between the plate and the 
pipe is neglected. The model can be extended by modeling the pipe as well 
similar to the other components of the collector. 



APPENDIX A 


The values of the parameters used are listed below : 


Dimension of the collector 

1.15 mxl.05 m 

Width of plate 

0.95 m 

Length of.Plate 

1 m 

Thickness of the glass cover 

4 mm 

Thickness of absorber plate 

1mm 

Thickness of insulation 

80 mm 

Inner diameter of pipe 

10.2 mm 

outer diameter of pipe 

12.5 mm 

Material properties of the collector 

Thermal conductivity of the glass covers 

.078 W m-i K-i 

Thermal conductivity of the copper absorber plate 

376. W m-^ K-i 

Thermal conductivity of the glass wool insulation 

.038 W m-^ K-^ 

Thermal conductivity of the fluid(water) 

.59 W m-i K-i 

Specific heat of the glass covers 

840 J kg-^ K-i 

Specific heat of the absorber plate 

383.1 J kg-^ K-i 

Specific heat of the insulation 

700 J kg-i K-^ 

Specific heat of the fluid 

4187 J kg-i K-i 

Density of the glass cover 

2700 kg m-3 
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Density of the absorber plate 

8954 kg m“® 

Density of the insulation 

24 kg 

Density of the fluid 

1000 kg 

Refractive index of glass relative to air 

1.526 

Extinction coefficient of glass 

6.8504 m-i 

Reflectivity of the surrounding 

0.2 

Latitude of kanpur 

26.28° 

Discretization grids 

Cover2 

11x11 

Coverl 

11x11 

Absorber plate 

11x11 

Insulation 

11x11x11 

Fluid 

11 

At 

10 s 
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